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ABSTRACT 

Models of structure formation with a cosmological constant A provide a good fit to the 
observed power spectrum of galaxy clustering. However, they suffer from several problems. 
Theoretically, it is difficult to understand why the cosmological constant is so small in Planck 
units. Observationally, while the power spectra of cold dark matter plus A models have ap- 
proximately the right shape, the COBE-normalized amplitude for a scale invariant spectrum 
is too high, requiring galaxies to be anti-biased relative to the mass distribution. Attempts 
to address the first problem have led to models in which a dynamical field supplies the vac- 
uum energy, which is thereby determined by fundamental physics scales. We explore the 
implications of such dynamical A models for the formation of large-scale structure. We find 
that there are dynamical models for which the amplitude of the COBE-normalized spectrum 
matches the observations. We also calculate the cosmic microwave background anisotropies 
in these models and show that the angular power spectra are distinguishable from those of 
standard cosmological constant models. 



1 Introduction 



The cosmological constant has had a long and tortured history since Einstein first intro- 
duced it in 1917 in order to obtain static cosmological solutions 0. Under observational 
duress, it has been periodically invoked by cosmologists and then quickly forgotten when 
the particular crisis passed. Historical examples include the first 'age crisis' arising from 
Hubble's large value for the expansion rate (1929), the apparent clustering of QSO's at a 
specific redshift (1967), and early cosmological tests which indicated a negative deceleration 
parameter (1974). 

Recently, a cosmological model with substantial vacuum energy — a relic cosmological 
constant A — has again come into vogue for several reasons @. First, dynamical estimates of 
the mass density on the scales of galaxy clusters, the largest gravitationally bound systems, 
suggest that Q m = 0.2 ± 0.1 for the matter (m) which clusters gravitationally (where Q is 
the present ratio of the mean mass density of the universe to the critical Einstein-de Sitter 
density, Q = 87tGp/3H 2 ) ||. However, if a sufficiently long epoch of inflation took place 
during the early universe, the present spatial curvature should be negligibly small, Qtot — 1- 
A cosmological constant, with effective density parameter Q\ = A/3Hq = 1 — Q m , is one 
way to resolve the discrepancy between Q m and Q t ot- 

The second motivation for the revival of the cosmological constant is the 'age crisis' 
for spatially flat Q m = 1 models. Current estimates of the Hubble expansion parame- 
ter from a variety of methods appear to be converging to H ~ 70 ± 10 km/sec/Mpc, 
while estimates of the age of the universe from globular clusters are holding at t gc ~ 
13 — 15 Gyr or more. Thus, observations imply a value for the the 'expansion age' Hob® = 
(H /70 km/sec/Mpc) (t /14 Gyr) ~ 1.0 ± 0.2. This is higher than that for the standard 
Einstein-de Sitter model with Q m = 1, for which H t = 2/3. On the other hand, for models 
with a cosmological constant, Hoto can be larger. For example, for = 0.6 = 1 — Q m , 
H t = 0.89. 

Third, cold dark matter (CDM) models for large-scale structure formation which include 
a cosmological constant (hereafter, ACDM) provide a better fit to the shape of the observed 
power spectrum of galaxy clustering than does the 'standard' fl m = 1 CDM model [|J. 
Figure 1 shows the inferred galaxy power spectrum today (based on a recent compilation 
||), compared with the matter power spectra predicted by standard CDM and a ACDM 
model with Q\ = 0.6. In both cases, the Hubble parameter has been fixed to h = ifo/(100 
km/sec/Mpc) = 0.7 and the baryon density to Qb = 0.0255, in the center of the range allowed 
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by primordial nucleosynthesis. Linear perturbation theory has been used to calculate the 
model power spectra, P(k), defined by (<5(k)5*(k')) = (27r) 3 P(fc)<5£)(k — k'), where <5(k) is 
the Fourier transform of the spatial matter density fluctuation field and 5d is the Dirac delta 
function. Here and throughout, we have taken the primordial power spectrum to be exactly 
scale- invariant, P pr imordiai{k) oc k n with n = 1. Standard CDM clearly gives a poor fit to the 
shape of the observed spectrum |J, while the ACDM model gives a good fit to the shape 
of the observed spectrum. The amplitudes of the model spectra in Fig. 1 have been fixed 
at large scales by observations of cosmic microwave background (CMB) anisotropies by the 
COBE satellite [g |. 

Despite these successes, cosmological constant models face several difficulties of their 
own. On aesthetic grounds, it is difficult to understand why the vacuum energy density of 
the universe, pa = A/8ttG, should be of order (10~ 3 eV) 4 , as it must be to have a cosmological 
impact (Q\ ~ 1). On dimensional grounds, one would expect it to be many orders of 
magnitude larger - of order mp lanck or perhaps mg USY - Since this is not the case, we might 
plausibly assume that some physical mechanism sets the ultimate vacuum energy to zero. 
Why then is it not zero today? 

The cosmological constant is also increasingly observationally challenged. Preliminary 
results from on-going searches|5| for distant Type la supernovae indicate that f^A < 0.47 (at 
95% confidence) for spatially flat A models. Furthermore, in A models a larger fraction of 
distant QSOs would be gravitationally lensed than in a A = universe; surveys for lensed 
QSOs have been used to infer the bound Q\ ~ 0.7 fllPf . 

In this paper, we focus on a third problem of cosmological constant models — the ampli- 
tude of the power spectrum of galaxy clustering. The shape of the ACDM power spectrum 
in Figure 1 matches the galaxy power spectrum; however the amplitude is too high. Indeed 
a number of analyses have found that this problem persists on all scales: 

• On the largest scales (k < O.lh Mpc -1 ), linear theory should be adequate, and Figure 
1 suggests that the amplitude is too high by at least a factor of two. 

• On intermediate scales, we can quantify the amplitude through the dispersion of the 
density field smoothed over top-hat spheres of radius R = 8/i _1 Mpc, denoted <r 8 , 
where a 2 (R) = 4-n f£° k 2 P(k)W 2 (kR)dk, and W(kR) is the Fourier-transform of the 
spatial top-hat window function of radius R. In the ACDM model of Figure 1, COBE 
normalization yields a 8 ~ 1.3 ||, while galaxy surveys generally indicate a Sjga i — 1 
for optically selected galaxies and ~ 0.8 for galaxies selected by infrared flux. This 
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high COBE normalization also marginally conflicts with the abundance of rich galaxy 



clusters JTTJ. Using the observed cluster X-ray temperature distribution function and 
modelling cluster formation using Press-Schechter theory, for this ACDM model the 
cluster abundance implies erg ~ l.Ojl;!! H3' wriere the errors are approximate 95% 
confidence limits. 



N-body simulations indicate that the power spectrum amplitude is higher by a factor 
of two to three than that found in galaxy surveys at small scales, k > OAh Mpc -1 \\L3\ . 
Thus, the cosmological constant model would require galaxies to be substantially anti- 
biased with respect to the mass distribution, a ga i < a p . Models of galaxy formation, 



however, suggest that the bias parameter, b = a ga i/a p , is greater than unity [[14], |i5 



Motivated by these difficulties, we consider models in which the energy density resides in a 
dynamical scalar field rather than in a pure vacuum state. These dynamical A models |16|, |I~7| 
were proposed in response to the aesthetic difficulties of cosmological constant models. They 
were also found |16| to partially alleviate their observational problems as well; for example, 
the statistics of gravitationally lensed QSOs yields a less restrictive upper bound on H t 



in these modelsf^]. We emphasize here that they may also solve the galaxy clustering 
amplitude problem. 

To get a preview of this conclusion, Fig. 1 also shows the COBE normalized power spec- 
trum for a dynamical A model with present scalar field density parameter = 0.6 (see §3 
for a discussion of these models). While the shape of the spectrum is identical to that of the 
ACDM model with Q\ = 0.6, the scalar field model has a lower amplitude, and thus provides 
a better fit to the galaxy clustering data. In §2, we explain these features of the power spec- 
trum for the standard ACDM model and for generic dynamical A models. The remaining 
sections investigate in detail a specific class of models as a worked example. Section 3 re- 
views the scalar field model, based on ultra-light pseudo-Nambu-Goldstone bosons (PNGBs) 
16fl . To explore the parameter space of this model, we have adapted a code which solves 



the linearized Einstein-Boltzmann equations for perturbations to a Friedmann-Robertson- 
Walker (FRW) background. The appendices contain details of these modifications. Section 
4 discusses the qualitative features of cosmic evolution in the PNGB models and presents 
results of our calculation for the amplitude of the power spectrum in this model. In §5 we 
present the cosmic microwave background (CMB) power spectrum for a particular set of 
model parameters, followed by the conclusion. 
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2 The Power Spectrum 

2.1 The Shape of P{k) 

Figure 1 suggests that standard CDM could be improved by simply shifting the turnover in 
the power spectrum to larger scales (smaller wavenumber k). This is a plausible fix, for the 
location of the turnover corresponds to the scale that entered the Hubble radius when the 
universe became matter-dominated. On scales smaller than this, the fluctuation amplitude is 
suppressed compared to that on larger scales, because matter perturbations inside the Hubble 
radius cannot grow in a radiation-dominated universe. This scale is determined by the ratio 
of matter to radiation energy density at early times. To "fix" CDM, one must decrease the 
ratio prn/pr in the universe today below that predicted by the standard Einstein-de Sitter 
model. The matter and radiation densities scale as p m = p m ,oa -3 and p r = p rj oa -4 , where the 
cosmic scale factor a is normalized to unity today (a = 1) and the subscript denotes the 
present. Thus the epoch of matter-radiation equality is determined by the present energy 
densities of matter and radiation: 

p r , 4.3 xlO- 5 
a EQ = — - Q p ■ W 

Decreasing the matter to radiation density ratio shifts the epoch of matter-radiation equality 
closer to the present, thereby moving the turnover in the power spectrum to larger scales. 

Indeed, this shift is precisely what is done in several currently popular models of structure 
formation. Examples include i) models with a lower Hubble constant than indicated by 
observations ||19||, ii) models with extra relativistic degrees of freedom|2(|, and iii) models 



with a cosmological constant H. Since p m oc Q m h 2 , a lower Hubble constant decreases the 
ratio of matter to radiation density today. Adding more relativistic degrees of freedom adds 
to the radiation content, decreasing the ratio of matter to radiation. Finally, in spatially 
flat A models, Q m = 1 — is reduced from its standard CDM value (0 OT = 1), achieving a 
similar effect. 

Thus the main benefit of A models for the shape of the power spectrum is that Q m is 
smaller than in the standard CDM model. For the purpose of the power spectrum shape, 
the value of the vacuum energy density at early times is irrelevant, as long as it is negligible 
compared to the matter and radiation densities at matter-radiation equality. While the time 
dependence of the vacuum energy density is different for various dynamical A models, all 
such models yield the same power spectrum shape for a fixed value of the present vacuum 
energy density. We emphasize this point in Fig. 2, which shows the energy densities of 
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matter, radiation, A, and a specific dynamical A model (scalar field <fi), as a function of scale 
factor a. With Q\ and = 0.6 today, the standard and dynamical A models have the same 
shape for P(k) (shown in Fig. 1), since they have identical values of oeq. As we will see 
below, however, the amplitudes of the power spectra in these models differ substantially. 

2.2 The Amplitude of P(k) 



Compared to standard CDM, three new physical effects |2T| conspire to change the ampli- 
tude of the matter power spectrum in COBE-normalized A models: i) the suppression of 
growth of perturbations when the universe becomes A-dominated, ii) the reduced gravita- 
tional potential, and iii) the integrated Sachs- Wolfe (ISW) effect. We review these effects in 
turn. 

The equations governing large scale perturbations in a flat universe with matter and 
vacuum energy are 

5 + Ha5 - ^H 2 Q m 5 = (2) 



H 2 - H ° 



a 3 



f2 m + Vl/y— — o? 



(3) 



Pa,o 

Here overdots denote derivatives with respect to conformal time r, where r = J dt/a(t), 
Pa is the vacuum energy density, not necessarily equal to its present value pa,o, the density 
fluctuation amplitude 5(x, r) = (p m (x, r) — p m (r)) / p m (r) , and H is the Hubble expansion 
rate [we use units in which H = c = 1]. 

Equation essentially describes the behavior of a damped harmonic oscillator. When the 
energy density of the universe becomes dominated by a A or dynamical A, i.e., the second 
term on the the RHS in equation becomes important, the damping becomes more severe. 
When this happens, the growth of perturbations is suppressed. As a function of Q m , this 
suppression can be described by the scaling, 

5 /5(z=wo) oc tt p m , (4) 

where 5o is the perturbation amplitude today, and 5( z =ioo) is the amplitude at the epoch 
z = (I/a) — 1 = 100, chosen as an arbitrary early epoch before the vacuum energy becomes 
dynamically important. In ACDM models, p ~ 0.2. For dynamical A models, the suppression 
exponent depends on the details of the specific model, but it is generally greater than that in 
ACDM models, because the dynamical A dominates earlier in the history of the universe for 
fixed Pa,o- F° r t ne model shown in Fig. 1, p — 0.56. For open CDM models (with A = 0), 
the scaling is also p ~ 0.56. 
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As a result of this suppression, one might expect the amplitude of the power spectrum in 
ACDM and dynamical A models to be smaller than that in standard CDM. However, from 
the Poisson equation 

V 2 $ = ^H 2 n rn 5 , (5) 

we have $ oc Q m 5, where $ is the gravitational potential associated with large-scale density 
fluctuations. Since the CMB anisotropy at large angle is a well-defined function of the 
potential^], COBE normalization corresponds to fixing the potential, i.e., to fixing Q m 5. 
For COBE-normalized models, the growth suppression and Poisson's equation combine to 
yield the scale-independent relation 5 oc fi^ 1 . Thus the power spectrum P(k) oc 5 2 oc fi" 1 ' 6 
in ACDM models. A larger cosmological constant implies a smaller Q m , which in turn implies 
a larger amplitude for the power spectrum. In dynamical A models, p is not fixed at 0.2, 
so the amplitude of the power spectrum can be smaller than in standard A models. For the 
model of Fig. 1, with p = 0.56, P(k) oc fi~°- 9 . 

The integrated Sachs- Wolfe effect (ISW), which is due to time evolution of the potential, 
also affects the amplitude of the power spectrum. The changing potential at late times in 
A models increases the anisotropy on the large angular scales probed by COBE. Thus, for 
fixed COBE normalization, the amplitude of the power spectrum decreases, changing the 
dependence of the power spectrum on Q m to P oc fi^ 1 ' 4 i n the ACDM model. In dynamical 
A models, where the potential typically changes more than in standard A models, the ISW 
effect tends to be larger and is not a power law function of Q m . Hence the power spectrum 
amplitude in dynamical models is even less enhanced than in ACDM models, and can even 
be reduced compared to standard CDM. 



3 Ultra-light Scalar Fields 



A number of models with a dynamical A have been discussed in the literature fP7 |. We will 
focus on a particular class of models motivated by the physics of pseudo-Nambu-Goldstone 
bosons (hereafter PNGBs) ||T6| , [23 . 



It is conventional to assume that the fundamental vacuum energy of the universe is zero, 
owing to some as yet not understood mechanism, and that this mechanism 'commutes' with 
other dynamical effects that lead to sources of energy density. This is required so that, e.g., 
at earlier epochs there can temporarily exist non-zero vacuum energy which allows inflation 
to take place. With these assumptions, the effective vacuum energy at any epoch will be 
dominated by the heaviest fields which have not yet relaxed to their vacuum state. At late 
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times, these fields must be very light. 

Vacuum energy is most simply stored in the potential energy V(<fi) ~ M 4 of a scalar 
field, where M sets the characteristic height of the potential, and we set V(4> m ) = at the 
minimum of the potential by the assumptions above. In order to generate a non-zero A at 
the present epoch, must initially be displaced from the minimum (0j 7^ m as an initial 
condition), and its kinetic energy must be small compared to its potential energy. This 



implies that the motion of the field is still overdamped, = J \V"((f)i)\ < 3H = 5 x 10 33 h 
eV. In addition, for ~ 1, the potential energy density should be of order the critical 
density, M 4 ~ 3iJ 2 M^/8vr, or M ~ 3 x l(T 3 /i 1/2 eV. Thus, the characteristic height and 
curvature of the potential are strongly constrained for a classical model of the cosmological 
constant. 

This argument raises an apparent difficulty for such a model: why is the mass scale 
thirty orders of magnitude smaller than Ml In quantum field theory, ultra-low-mass scalars 
are not generically natural: radiative corrections generate large mass renormalizations at 
each order of perturbation theory. To incorporate ultra-light scalars into particle physics, 
their small masses should be at least 'technically' natural, that is, protected by symmetries, 
such that when the small masses are set to zero, they cannot be generated in any order of 
perturbation theory, owing to the restrictive symmetry. 

From the viewpoint of quantum field theory, PNGBs are the simplest way to have natu- 
rally ultra-low mass, spin-0 particles. PNGB models are characterized by two mass scales, 
a spontaneous symmetry breaking scale / (at which the effective Lagrangian still retains the 
symmetry) and an explicit breaking scale /i (at which the effective Lagrangian contains the 
explicit symmetry breaking term). In terms of the mass scales introduced above, generally 
M ~ /1 and the PNGB mass ~ /x 2 / /. Thus, the two dynamical conditions on and M 
above essentially fix these two mass scales to be /1 ~ M ~ 10 -3 eV, interestingly close to the 
neutrino mass scale for the MSW solution to the solar neutrino problem, and / ~ Mp\ ~ 10 19 
GeV, the Planck scale. Since these scales have a plausible origin in particle physics models, 
we may have an explanation for the 'coincidence' that the vacuum energy is dynamically 
important at the present epoch. Moreover, the small mass is technically natural. 

An example of this phenomenon is the 'schizon' model | 23| , based on a Z/v-invariant 
low-energy effective chiral Lagrangian for N fermions, e.g., neutrinos, with mass of order M, 
in which the small PNGB mass, ~ M 2 /f, is protected by fermionic chiral symmetries. 
The potential for the light scalar field is of the form 

y(0) = M 4 [cos(0//) + l] . (6) 
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Since is extremely light, we assume that it is the only classical field which has not yet 
reached its vacuum expectation value. The constant term in the PNGB potential has been 
chosen to ensure that the vacuum energy vanishes at the minimum of the potential, in 
accord with our assumption that the fundamental vacuum energy is zero. 

4 Cosmic Evolution and Large-scale Power Spectrum 
in PNGB Models 

To study the cosmic evolution of these models, we focus on the spatially homogeneous, 
zero-momentum mode of the field, <f>(°'(r) = (0(x, r)) 5 where the brackets denote spatial 
averaging. We are assuming that the spatial fluctuation amplitude #0(x, r) is small compared 
to as would be expected after inflation if the post-inflation reheat temperature Tjuj < 
f ~ M p i. The scalar equation of motion is given in Appendix A. 

The cosmic evolution of is determined by the ratio of its mass, ~ M 2 /f, to the 
instantaneous expansion rate, H(t). For < 3H, the field evolution is overdamped by 
the expansion, and the field is effectively frozen to its initial value 0j. Since is initially 
laid down in the early universe (at a temperature T ~ / 3> M) when its potential was 
dynamically irrelevant, its initial value in a given Hubble volume will generally be displaced 
from its vacuum expectation value m = nf (vacuum misalignment). Thus, at early times, 
the field acts as an effective cosmo logical constant, with vacuum energy density and pressure 

~ — ~ M 4 . At late times, 3> 3H(t), the field undergoes damped oscillations 
about the potential minimum; at sufficiently late times, these oscillations are approximately 
harmonic, and the stress-energy tensor of averaged over an oscillation period is that of 
non-relativistic matter, with energy density ~ a~ 3 and pressure ~ 0. 

Let t x denote the epoch when the field becomes dynamical, = 3H(t x ), with corre- 
sponding redshift 1 + z x = l/a(r x )) = (M 2 /3H f) 2 / 3 . For comparison, the universe makes 
the transition from radiation- to matter-domination at z eq ~ 2.3 x 10 4 Q m h 2 , much earlier 
than when the field becomes dynamical. The f — M parameter space is shown in Fig. 3. 
In the far right portion of the figure, the field becomes dynamical before the present epoch 
and currently redshifts as non-relativistic matter; on the far left, is still frozen and acts 
as an ordinary cosmological constant. In the dynamical region, the present density param- 
eter for the scalar field is approximately ~ 2An(f /Mpi) 2 , independent of M ||24j| . The 



quasi- horizontal lines show contours of constant fl^, assuming a typical initial field value 
0i// = 1.6 (we will use this value of (pi/ f for all the plots below; the quoted limits and 
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results depend slightly on it). The limit < 1 corresponds approximately to / < 3.5 x 10 18 
GeV. In the frozen region, on the other hand, is determined by M 4 , independent of /, and 
the contours of constant Q<p are nearly vertical. In this region, the bound < 1 corresponds 
roughly to M < 0.003 eV. 

Figure 4 shows contours of constant H t i n the same parameter space. As expected, 
models with large ifo^o are concentrated toward the left hand portion of the figure; as one 
moves to the right, H t asymptotically approaches the Einstein-de Sitter value 2/3, since 
the scalar field currently redshifts as non-relativistic matter and we have assumed a spatially 
flat universe. Consequently, the 'interesting' region of parameter space is the area near the 
'corner' in Figs. 3 and 4, in which the field becomes dynamical at recent epochs, z x ~ — 3. 
This has new consequences, compared to A models, for the classical cosmological tests, the 
expansion age H t , and large-scale structure. In this region, the mass of the PNGB field is 
miniscule, ~ 3i?o ~ 4 x 10~ 33 eV, and (by construction) its Compton wavelength is of 
order the current Hubble radius, = m^ 1 = Hq 1 /3 ~ 1000/iT 1 Mpc. 

Figure 5 shows contours of the amplitude of galaxy clustering in the f — M parameter 
space. The amplitude shown is the quantity 

~{P{k)/k)i 



lim 



(7) 



(P(k)/k) A _ 

i.e., the amplitude on large scales relative to that for a ACDM model with the same effective 
density as the PNGB model, f^A = This amplitude ratio goes to unity in the left-hand 
portion of the figure since that region corresponds to a ACDM model. However the amplitude 
ratio can be substantially below one in the dynamical region on the right. The cross marks 
the specific choice M = 0.005 eV, / = 1.885 x 10 18 GeV, with initial field value fc/f = 1.6, 
yielding Vt^ = 0.6, which corresponds to the parameters used for the dynamical A curves 
in Figs. 1 and 2. For this case, the X-ray cluster abundance yields of ~ 0.9j/|, in good 
agreement with the COBE normalization a% OBE ~ 0.8 for this model. Figure 6 shows how 
density perturbations grow in the different models. From Eqn.(4) and the text following, 
the dynamical A model has a higher amplitude at early times than a ACDM model with 
the same amplitude today. As a consequence, there should be no problem accounting for 
high-redshift objects such as QSOs and Lyman-alpha clouds in this model. 

Note that the factor 5(z)/5o, relative to its value in the standard CDM model, approaches 
0~ p at z ^> 1, where p is the scaling exponent discussed in §2. As a result, the non-linear 
behavior of the dynamical A model follows that of an open model with the same value of 



fl m . We estimate the non-linear behavior by using the fitting formula of Ref. [25], following 
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the original treatment |26| of Hamilton et al. Figure 7 shows these non-linear spectra. On 
scales k < 1/iMpc -1 , the amplitude of the power spectrum is indeed a factor of two smaller 
in the dynamical A model than in the corresponding ACDM model. 

We note by comparing Figs. 4 and 5 that the region of parameter space in which the 
amplitude (anti-bias) problem is solved, i.e., in which the amplitude ratio is approximately 
in the range 0.3 — 0.5, is the one in which the age of the universe is only slightly greater 
than in the Einstein-de Sitter Q m = 1 case. For our specific model above, -f^o = 0.73. For 
the corresponding ACDM model with the same value of Q m , H to = 0.89, more comfortably 
within the observational limits. This is a general feature of the dynamical models considered 
here: for fixed fl m , the standard A model gives an upper bound on H t . Thus, the amplitude 
problem in this model is resolved partially at the expense of the age problem. On the other 
hand, the go constraints from SNe and gravitational lensing translate into weaker upper 
bounds on H t for the dynamical as opposed to the standard A models. Although we 
have not thorougly examined all models, it is clear that one could explore the PNGB model 
parameter space to obtain a more balanced compromise between the age problem and the 
anti-bias problem. For example, for / ~ 2.5 x 10 18 GeV and M ~ 0.0035 eV, the amplitude 
ratio is about 0.5, and one has ~ 0.75 and H^q ~ 0.9. In this case, with h = 0.7, the 
power spectrum shape is reasonable (Sl m h ~ 0.15) and the age of the universe is t ~ 12.6 
Gyr. 

Comparing Figs. 3 and 5, and focusing on the dynamical region near the 'corner' of the 
parameter space, we see that the power spectrum shape and amplitude constraints fix the 
free parameters of the model. That is, as noted in §2, the shape of the spectrum is fixed by 
requiring VL^ ~ 0.6, which determines the scale /. Near the corner, fixing the amplitude then 
determines the other mass scale M. While these figures correspond to a specific choice of the 
initial field value (f>i/f, the scalar field evolution is universal in the sense that a shift in the 
mass scale /, accompanied by an appropriate rescaling of 0j, leads to essentially identical 
evolution. Consequently, compared to ACDM models, these dynamical models have only 
one additional free parameter, the mass M, to solve the amplitude (anti-bias) problem. 

5 CMB Anisotropy 

The angular power spectra of the cosmic microwave background (CMB) anisotropy for dy- 
namical A models are distinguishable from those of standard CDM and ACDM models. 
CMB angular power is usually expressed in terms of the angular multipoles G\. If the sky 
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temperature is expanded in terms of spherical harmonics as T(9, 0) = E/ m az m Yj m (0, 0), then 
Ci = (\ai m \ 2 ), where large / corresponds to small angular scales. The angular power spec- 
tra for standard CDM (Q m = 1), ACDM, and dynamical A models (the latter two with 
VL, m = 0.4) are shown in Figures 8 and 9 for h = 0.7, £l B = 0.0255, and primordial spectral 
index n — 1. Following standard practice, we plot the product 1(1 + l)Cj, normalized to its 
value at / = 10, vs. I. 

The Appendices contain the details of the alterations required in the standard Boltzmann 
code to calculate the CMB anisotropy in scalar field dynamical A models. We can, however, 
identify two physical effects primarily responsible for the differences in the CMB signature 
between the ACDM and dynamical A models shown in Figs. 8 and 9. First, the present 
ages in conformal time coordinates, r , are different in the two models. Even though the 
acoustic oscillations responsible for the peaks in the CMB angular spectrum occur at the 
same physical scales (or same Fourier wave numbers k), the correspondence between k and 
angular multipole / differs. Typically, in a flat universe, a given multipole / corresponds to 
a fixed value of kr . Thus, the dynamical A angular spectra are shifted in I by the ratio 
of the present conformal times in the two models. Second, since the scalar field evolves at 
late times, the gravitational potential changes more rapidly in the dynamical A model. This 
leads to an enhanced ISW effect, and therefore a relatively larger C x at large scales (low /), 
as shown in Fig. 9. Thus, for models normalized by COBE, which approximately fixes the 
spectrum at I ~ 10, the angular amplitude 1(1 + l)Ci at small scales (large I) is smaller in 
the dynamical A model. 

6 Conclusions 

The observational arguments in favor of the resurrection of the cosmological constant apply 
to dynamical A models as well. In addition, the dynamical A models offer a potential 
physical explanation for the curious coincidence that Q\ is close to one, by relating the 
present vacuum energy density to mass scales in particle physics. In the ultra-light pseudo- 
Nambu-Goldstone boson models, this is achieved through spontaneous symmetry breaking 
near the Planck scale, / ~ M P i, and explicit breaking at a scale reminiscent of MSW neutrino 
masses, M ~ 10 -3 eV. In combination with the assumption that the true vacuum energy 
vanishes (due to an as yet unknown physical mechanism), such a model provides an example 
of a dynamical A. 

We have shown that such dynamical models can lead to a lower amplitude for density 
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fluctuations compared to standard A models, thereby alleviating the anti-bias problem. The 
advantages of the cosmological constant for the shape of the power spectrum are retained in 
the dynamical models as well. Such dynamical models are, moreover, distinguishable from 
constant-A models by virtue of their CMB angular spectra. 
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A Changes in standard Boltzmann code 

This Appendix and the following briefly outline the new physics incorporated into the Boltz- 
mann code in ACDM and dynamical A models. Since the Hubble parameter is determined 
by the sum over densities of all species, H 2 = (8nG)T>iPi, inclusion of a cosmological constant 
A or scalar field <fi changes the relationship between the cosmic scale factor a and conformal 
time r, since (da/dr)/a 2 = H. In addition to the species included in the standard Boltzmann 
code, namely, baryons, cold dark matter, photons, and three massless neutrinos, the density 
in a cosmological constant or scalar field <fi is now included. In ACDM models, the vacuum 
energy density p\ = A/87rG is constant. In the dynamical models, the scalar field energy 
density p^ can be solved for with the scalar equation of motion for the homogeneous part 
0(°)(r) of the field, 

+ 2Haj>W + a?dV(<f><®)/d(/><® = , (8) 
where the scalar field potential is 

V(0) = M 4 [cos(0//) + l] , (9) 

and the scalar energy density 

P , = ^-J m + v(^) . (10) 

Here overdots denote derivatives with respect to conformal time r. 
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B Perturbation equations for dynamical A models 



The general equation of motion for the scalar field 0(x, r) is derived by minimizing the action 



S 



-g^d^-V^) 



(11) 



with respect to variations in 0. The metric is that of a perturbed Friedmann-Robertson- 
Walker universe, 

^(x,r) =^(r) + %,(x,r) , (12) 

where g^J is the homogeneous part which describes the Hubble expansion, and 5g^ v is the 
metric perturbation. In synchronous gauge, the latter can be parametrized by the variables 
h, hz3 as in [pT]] . The scalar field can be similarly decomposed into a homogeneous part and 
a spatial perturbation, 



0(x,r) = 0(°)(r) + ( 50(x,r) , (13) 

where is the solution to the spatially homogeneous equation of Appendix A. Keeping 
only terms linear in h, and 5<fi, and taking the Fourier transform yields the equation of 
motion for the Fourier amplitude 5<pk, 



+ 2Ha(50 k ) + (k 2 + a 2 [d 2 V/d(f) 2 



I 0=0(0) ( T ) 



hj><® 



(14) 



There will also be an additional source term in the Einstein equation for the metric 



perturbation. Again following the notation of [27], the Einstein equation becomes: 



h + Hah = 8ttG (£0 + S u ) 
where the source term due to <fi is given by 

S+ = 4(50)0 - 2a 2 (5<P)[dV/d<P} (f)=4>{0){T) , 



(15) 



(16) 



and S n , contains the usual source terms for matter and radiation [27 . 
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Figure V. COBE-normalized power spectra for standard CDM, ACDM with VL^ = 0.6, and 
scalar field = 0.6 models. In all models h = 0.7, Qb — 0.0255, and n = 1. The data 
points are based on a recent compilation of galaxy clustering data by Peacock and Dodds 
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Figure 2: Density p vs. cosmic scale factor a. Fixing I^a or to 0.6 lowers Vt m from the 
standard CDM value of 1.0, pushing the epoch of matter-radiation equality, oeq, closer to 
today. The cross denotes a E Q for the standard CDM model and the asterisk denotes a E Q 
for the A and models. 
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Figure 3: Contours of in the PNGB parameter space, assuming an initial field value 
<pi/f = 1.6. The cross marks the choice M = 0.005 eV, / = 1.885 x 10 18 GeV, yielding 
Qcf, = 0.6, which is the model shown in Figures 1, 2, and 6-9. 
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Figure 4: Contours of H t in the PNGB parameter space for (pi/ f = 1.6. The cross indicates 
the same model as in Fig. 3. 
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Figure 5: Amplitude contours in the PNGB parameter space for fa/ f = 1.6. The amplitude 
shown is defined as lim^ 



(P(k)/k)_ 



, the amplitude on large scales relative to that of a 
ACDM model with the same effective density as the PNGB model, Qa = Q<f>- Again, the 
cross marks the sample model for which both the power spectrum shape and amplitude 
provide a good fit to the galaxy clustering data. 
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Figure 6: Evolution of density perturbations. Shown is the density fluctuation amplitude 
at redshift z normalized to its present amplitude, S(z)/5q, vs. z. The models shown are 
standard Q m = 1 CDM (solid-red), ACDM with Q A = 0.6 = 1 - Q m (dotted-green), an 
open CDM model with Q m = 0.4 (short dashed-light blue), and the dynamical A model with 
flff, = 0.6 (long dashed-dark blue ). 



22 




Figure 7: Power spectra for COBE- normalized standard CDM (solid-red) with cr 8 = 1.2; 
ACDM with Q\ = 0.6 (dashed-blue), for which <7g = 1.0; and the dynamical A model with 
flff, = 0.6 (dotted-green), for which a 8 = 0.8. The latter two models are normalied to the 
cluster abundance and have h = 0.7. Lower curves show the linear theory power spectra, 
upper curves the non-linear spectra obtained from scaling relations extracted from N-body 
simulations. 
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Figure 8: CMBR angular power spectra for standard CDM, ACDM with — 0.6, and 
scalar field = 0.6 models. In all models h = 0.7, Qb = 0.0255, and n — 1. Plotted is 
1(1 + l)Ci vs. /, normalized at I = 10. 
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Figure 9. 



CDM 
A =0.6 
0» = 0.6 




I , , , I , , , I , , , I , , , I , , , I 

20 40 60 80 100 



Figure 9: Same as Fig. 8, but showing only the low I multipoles to emphasize the enhanced 
ISW effect in the PNGB model. 
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